function [meanerr,varerr]=checkou(alpha,m,y0,beta,pathnum)

t=1;
mcmean=0; mcvar=0;
for ii=1:pathnum
    currentpath=ouprocess(alpha,m,y0,beta,10000);
    mcmean=mcmean+currentpath(10000);
    mcvar =mcvar +currentpath(10000)^2;
end
mcmean=mcmean/pathnum; mcvar=mcvar/pathnum-mcmean^2;
exmean=y0*exp(-alpha*t)+m*(1-exp(-alpha*t));
if alpha ~=0
exvar =beta^2/2/alpha*(1-exp(-2*alpha*t));
else
exvar =1;
end
meanerr=abs(exmean-mcmean);
varerr =abs(exvar -mcvar );

